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ABSTRACT 


An  experimental  investigation  of  the  simulation  of  the  environment 
encountered  in  the  lower  regions  of  the  ionosphere  (50  to  500  km)  has 
been  conducted  at  AEDC.  Microwave  cavities  were  utilized  in  conjunc¬ 
tion  with  a  converging-diverging  nozzle  and  a  cryogenically  pumped 
vacuum  chamber  to  produce  supersonic  flow  fields  with  small  ion  num¬ 
ber  densities.  A  platinum,  cylindrical,  electrostatic  probe  was  used  to 
measure  the  voltage -current  characteristics  of  the  low  density  plasma 
in  the  nozzle  exit  plane.  For  exit  plane  Mach  numbers  of  2.86  to  3.  17, 
ion  densities  of  5.  75  x  10®  cm-^  to  1.  69  x  10®  cm“®  and  electron  tem¬ 
peratures  of  3000°K  ±900°K  were  calculated  using  the  electrostatic 
probe  characteristics. 
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SECTION  I 
INTRODUCTION 


In  recent  years  there  has  been  considerable  interest  (for  example, 
from  the  Air  Force  Cambridge  Research  Laboratory)  in  a  ground-based 
facility  capable  of  generating  conditions  corresponding  to  those  in  the 
ionosphere.  Such  a  facility  could  be  conveniently  utilized  in  the  calibra¬ 
tion  of  atmospheric  electrostatic  sampling  probes  which  attempt  to  deter¬ 
mine  ion  concentrations  and  electron  temperatures  in  the  upper  atmos¬ 
phere.  The  conditions  for  which  simulation  is  desired  are  (1)  altitudes 
from  50  to  500  km,  (2)  Mach  numbers  from  2  to  5,  and  (3)  ion  densities 
from  103  to  103  per  cc.  The  attempt  to  attain  these  conditions  simul¬ 
taneously  in  a  ground  test  facility  is  the  subject  of  this  report. 

The  upper  atmosphere  is  an  extremely  rarefied,  partially  ionized 
gas.  The  concentration  of  neutral  particles,  molecules,  and  atoms 
decreases  sharply  with  increasing  altitude  (Ref.  1).  For  example,  at 
100  km  there  are  approximately  10*3  particles/ cm3,  whereas  at  an  alti¬ 
tude  of  1000  km  there  are  approximately  103  particles /cm3  (Table  I, 
Appendix  n).  In  the  lower  regions  of  the  ionosphere,  the  concentrations 
of  ions  and  electrons  vary  strongly,  not  only  with  altitude  (Fig.  1,  Appen¬ 
dix  I)  but  also  with  the  time  of  year  and  during  a  24-hr  period  (Ref.  2).  . 
With  increasing  height,  the  concentrations  of  ions  and  electrons  decrease 
considerably  slower  than  the  molecular  concentration.  Accordingly,  the 
degree  of  ionization  of  the  ionosphere  increases  with  increasing  height. 
For  example,  at  an  altitude  of  100  km,  the  plasma  is  very  weakly 
ionized;  the  ratio  of  ions  to  neutrals  (N0/n0)  amounts  to  only  10  "8  to 
10 "10.  At  an  altitude  of  300  km  the  ratio  No/no*  10,  and  at  3000  km 
the  ratio  N0/hq  is  approximately  10^.  In  this  report  the  simulation  of 
the  lower  ionosphere  (50  to  500  km),  corresponding  to  conditions  which 
would  be  encountered  by  gages,  electrostatic  probes,  and  other  devices 
mounted  in  sounding  rockets,  is  discussed. 


SECTION  II 

EXPERIMENTAL  APPARATUS 


2.1  RESEARCH  VACUUM  CHAMBER  AND  CONICAL  NOZZLE 

The  vacuum  chamber  used  in  these  experiments  is  a  stainless  steel 
cylinder  4  ft  in  diameter  and  10  ft  long  (Fig.  2).  The  chamber  is  equipped 
with  a  300 -cfm  Stokes  mechanical  pump  which  is  used  to  rough  the  cell 
down  to  a  pressure  of  approximately  20  juHg.  Chamber  pressures  in  the 
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10“  torr  range  are  achieved  by  a  6 -in.  oil  diffusion  pump  backed  by  a 
13-cfm  Welch  mechanical  pump.  The  chamber  is  also  equipped  with  a 
liquid-nitrogen-cooled  cryoliner  (77°K)  and  a  finned  aluminum  gaseous- 
helium-cooled  cryopump  (20°K).  When  cryogenic  cooling  is  used,  a 
base  pressure  of  5  x  10"  ^  torr  is  achieved  in  the  chamber.  A  minia¬ 
ture  GE  ionization  gage  is  located  in  the  aft  portion  of  the  chamber  and 
is  used  to  monitor  the  static  pressure  in  this  region  (Fig.  2).  An 
Alphatron®  gage  is  located  in  the  plenum  chamber  and  is  used  to  record 
the  stagnation  pressure. 

A  stainless  steel  conical  Mach  3  nozzle  is  located  at  the  forward 
end  of  the  chamber  with  an  exit  diameter  of  20  in.  and  throat  diameter 
of  7.  1  in.  The  overall  length  of  the  nozzle  is  37  in.  (Fig.  3).  The 
throat  and  divergent  portion  of  the  nozzle  are  cooled  by  a  liquid- 
nitrogen  jacket. 


22  ION  GENERATION 

Three  Raytheon  PGM- 10  microwave  power  generators  featuring 
magnetron  oscillators  were  used  to  ionize  a  portion  of  the  gas  flow  into 
the  plenum  chamber.  Nitrogen  gas  from  a  high-pressure  bottle  was 
admitted  to  the  plenum  chamber  through  a  manifold  of  three  0.  5-in.-diam 
quartz  tubes  (Fig.  4).  Each  tube  passed  through  a  tunable  microwave 
cavity  where  energy  from  the  generators  was  applied  to  the  gas.  Each 
generator  is  capable  of  supplying  power  output  in  excess  of  85  watts  at 
a  frequency  of  2450  mHz.  Power  output  is  continuously  adjustable  by 
means  of  a  variable  auto -transformer  which  controls  the  voltage 
applied  to  the  magnetron. 


2.3  ELECTROSTATIC  PROBE 

The  electrostatic  probe  used  to  determine  ion  densities  was  a 
0.  508-cm-diam  platinum  rod  which  was  inserted  in  a  glass  tube,  leaving 
approximately  7  cm  of  its  length  exposed.  The  glass  tube  probe  support 
was  attached  to  the  side  of  the  pitot  probe  traversing  mechanism  within 
the  research  vacuum  chamber.  Teflon®  insulated  wire  was  used  to 
connect  the  probe  to  instrumentation  outside  the  chamber. 

A  block  diagram  of  the  electrostatic  probe  instrumentation  is  shown 
as  Fig.  5,  and  a  photograph  of  the  vacuum  chamber  and  instrumentation 
is  presented  in  Fig.  6.  Sweep  voltage,  furnished  by  a  Wavetec®  function 
generator,  was  a  triangular  wave  of  approximately  ±15  volts.  A  Boonton 
Model  95 A  sensitive  d-c  meter  was  used  to  measure  and  amplify  the 
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small  probe  currents.  The  meter  output,  a  voltage  proportional  to  input 
current,  was  applied  to  the  Y-axis  of  a  Moseley  X-Y  plotter.  Sweep 
voltage  from  the  function  generator  was  applied  to  the  X-axis,  producing 
current  versus  voltage  graphs  of  the  probe  characteristic  on  the  plotter. 
Isolation  plugs  were  used  with  all  instruments  so  that  the  only  ground 
connection  was  to  the  test  chamber. 


SECTION  III 

ELECTROSTATIC  PROBE  THEORY 


A  great  deal  of  information  regarding  the  nature  of  a  plasma  flow 
field  can  be  obtained  from  the  voltage- current  characteristics  of  an 
electrostatic  probe  immersed  in  the  charged  region.  The  following 
discussion  is  a  summary  of  the  probe  theory' presented  by  Langmuir 
and  Mott-Smith  (Ref.  3).  When  a  collector  is  negatively  charged,  it 
repels  electrons  from  its  neighborhood  but  gathers  positive  ions.  It 
thus  becomes  surrounded  by  a  positive  ion  sheath  which  contains  a  posi¬ 
tive  ion  space  charge  but  no  free  electrons.  The  whole  drop  in  potential 
between  the  ionized  gas  and  the  collector  becomes  concentrated  within 
this  sheath,  the  positive  space  charge  on  the  ions  in  the  sheath  being 
able  to  neutralize  the  effect  of  the  negative  charge  on  the  electrode  so 
that  the  field  of  the  collector  does  not  extend  beyond  the  outer  edge  of 
the  sheath.  The  number  of  ions  taken  up  by  the  collector  is  thus  limited 
by  the  number  that  reaches  the  outer  edge  of  the  sheath  as  a  result  of 
their  proper  motions.  The  current  density  over  the  area  of  the  outside 
of  the  sheath  measures  the  positive  ion  current  density  in  the  ionized 
gas. 


In  a  similar  manner,  a  positively  charged  probe  becomes  sur¬ 
rounded  by  an  electron  sheath,  and  the  current  is  limited  by  the  rate 
at  which  the  electrons  reach  the  edge  of  the  sheath;  and  in  this  manner, 
the  electron  current  density  in  the  ionized  gas  may  be  measured.  Be¬ 
cause  of  the  small  mass  and  high  mobility  of  electrons,  the  electron 
current  densities  in  uniformly  ionized  plasmas  are  hundreds  of  times 
greater  than  the  positive  ion  current  densities.  In  measuring  electron 
current  densities  there  is  a  limit  to  the  positive  voltage  that  may  be 
used  on  the  collector,  since  ionization  by  collision  becomes  significant 
at  a  certain  voltage.  Too,  the  positive  ions  produced  within  the  sheath 
neutralize  the  electron  space  charge  and  allow  the  current  to  increase 
indefinitely  without  further  increase  in  voltage. 

If  the  collector  is  at  a  potential  slightly  negative  with  respect  to  the 
space,  electrons  may  still  reach  the  probe  if  they  have  sufficient 
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velocities  to  carry  them  against  the  retarding  field  which  they  encounter 
within  the  positive  ion  sheath.  As  the  collector  is  made  more  negative, 
the  lower  speed  electrons  fail  to  reach  the  collector,  while  those  with 
high  speed  may  still  reach  it.  The  current-voltage  characteristic  of  a 
probe  thus  gives  indications  as  to  the  distribution  of  velocities  among 
the  electrons  in  the  ionized  gas. 

When  the  electron  velocities  have  a  Maxwellian  distribution,  the 
electron  current,  i,  flowing  to  any  collector  at  a  potential  V,  which  is 
negative  with  respect  to  the  surrounding  space,  is  given  by  the  equation 

i  =  IA  exp  {Ve/<kT>}  (1) 

where  I  is  the  electron  current  density  in  the  ionized  gas,  A  is  the  sur¬ 
face  area  of  the  probe,  e  is  the  charge  on  the  electron,  k  is  the 
Boltzmann  constant  (1.372  x  10“^-®  erg/°K),  and  T  is  the  absolute  tem¬ 
perature  corresponding  to  the  velocity  distribution  of  the  electrons.  If 
V  is  expressed  in  volts,  the  value  of  e/k  becomes  11,  600°K/volt. 

Equation  (1)  does  not  apply  when  the  potential  of  the  collector  is  such 
as  to  exert  an  attractive  force  on  the  particles  being  collected.  For  an 
accelerating  field,  the  apparent  increase  in  current  density  near  the 
probe  (in  the  sheath  region)  as  compared  with  that  farther  from  it,  is 
attributed  to  the  large  number  of  low -velocity  electrons,  resulting  from 
collisions,  which  are  trapped  within  the  region  around  the  collector  by 
the  accelerating  field  which  prevents  their  escape.  With  low  intensities 
of  ionization,  the  diameter  of  the  sheath  may  be  many  times  greater  than 
that  of  the  collector.  The  initial  velocities  of  the  electrons  or  ions  as 
they  enter  the  sheath  may  thus  cause  them  to  form  orbits  which  carry 
them  out  of  the  sheath  again  without  colliding  with  the  collector. 

When  the  sheath  radius  is  large  compared  to  that  of  the  collector, 
and  especially  when  the  initial  velocities  of  the  electrons  or  ions  are 
high,  the  orbital  motion  of  the  particles  must  be  considered.  Let  the 
small  circle  about  F  in  Fig.  7  represent  the  cross  section  of  a  cylindri¬ 
cal  probe  charged,  for  example,  negatively  with  respect  to  the  surround¬ 
ing  space.  The  field  produced  by  the  probe  extends  only  within  a  region 
of  radius  a,  indicated  in  the  illustration  by  the  dotted  circle.  Consider 
an  ion  of  mass  m  moving  with  a  velocity  v0  along  a  path  AB  and  let  BC 
be  the  extension  of  the  straight  line  AB.  Then  the  angular  momentum  of 
the  ion  about  F  is  equal  to  mvop0  where  pQ  is  the  perpendicular  distance 
OF  from  the  line  AC  to  the  axis  or  center  of  the  probe. 

The  force  exerted  on  the  ion  by  the  probe  will  cause  the  ion  to  move 
along  some  such  path  as  ABMDE,  the  portions  AB  and  DE  being  straight. 
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At  some  point  M,  the  ion  will  be  at  a  minimum  distance  pm  from  the 
center  of  the  probe,  and  at  this  point  the  velocity  vm  of  the  ion  is  wholly 
tangential.  Since  the  probe  can  exert  only  a  central  force  on  the  ion, 
the  angular  momentum  about  F  must  remain  constant  throughout  the  path. 
Thus,  one  may  place 

mvoPo  =  mvmpm  (2) 

It  is  convenient  to  express  the  velocities  of  ions  in  terms  of  the  potential 
differences  necessary  to  produce  them.  The  work  done  on  an  ion  in 
moving  through  a  potential  difference  V  is  Ve,  and  this  is  equal  to  the 
gain  in  kinetic  energy  so  that 

1/2  mv^  =  Ve  (3) 

Combining  this  with  Eq.  (2),  one  obtains 

Vo  P„2  -  VmPm2  (4) 

where  VQ  measures  the  initial  kinetic  energy  of  the  ion  in  volts  and  Vm 
is  the  corresponding  energy  when  at  the  minimum  distance  from  F.  As 
the  potential  is  increased  on  the  probe,  the  ion  is  deflected  further  so 
that  pm  decreases. 

When  pm  =  r,  the  radius  of  F,  the  ion  makes  a  grazing  collision 
with  the  surface  of  the  probe.  If  V  is  the  potential  of  the  probe,  then 
the  energy  of  the  ion  at  the  moment  of  collision  is  measured  by  VQ  +  V. 
Substituting  this  for  Vm  in  Eq.  4  and  placing  pm  =  r,  one  obtains 

Pn  =  r  VI  +  V/V„  (5) 

In  order  that  any  ion  may  reach  the  probe,  the  value  of  pQ  for  its 
path  must  be  less  than  the  value  of  Pomax  given  by  Eq.  5.  Thus,  p  is 

the  radius  of  the  ’’effective  target"  of  the  probe.  Therefore,  if  I  is  the 
current  density  of  the  ions,  the  current,  i,  taken  up  by  a  probe  of  sur¬ 
face  area  A  for  a  cylindrical  collector  is 

i  =  AI  y/1  +  V/VQ  (6) 

The  equations  are  equally  applicable  to  the  collection  of  electrons.  The 
problem  of  calculating  the  current  under  practical  conditions  is  more 
difficult  since  the  ions  do  not  all  have  the  same  initial  velocity  as  was 
assumed  in  deriving  Eq.  (6),  but  rather  have  velocities  distributed  by 
Maxwell's  law.  The  complete  mathematical  solution  of  this  problem  has 
been  worked  out  for  cylinders  by  Langmuir  and  Mott-Smith. 
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For  accelerating  fields,  the  solution  for  cylinders  of  surface  area 

A  is 


i  =  Alf 


(7) 


the  function  f  being  given  by 

f  =  (a/r)P  (yft)  +  er1{  1  -  P(y[W+^ )} 


where 


rj  =  Ve/{kT) 


2 


a  -  r 
2 

”  +  *  =~r—2 n 

a  -  r 

Here  P  stands  for  the  probability  integral  defined  by 

- 

The  symbols  V,  e,  k,  T,  a,  and  r  have  the  meanings  already  assigned. 

With  low  current  densities,  probes  of  small  radius,  and  high  initial 
velocities  of  the  ions,  the  sheath  radius  becomes  large  compared  to  that 
of  the  probe.  For  this  case,  one  may  obtain  from  Eq.  (7)  a  more  con¬ 
venient  equation  by  assuming  a/r  is  large  compared  to  unity  and  expand¬ 
ing  the  resulting  equation  as  a  reciprocal  power  series  in  a/r  which 
gives 


f  =  (2l^)yPT+l 


(8) 


Substituting  this  into  Eq.  (7)  and  squaring,  one  obtains 


.2 


l 


4  A2  I2 

7 T 


(Ve/(kT)  +  1) 


(9) 


Therefore,  if  the  square  of  the  observed  current,  i,  is  plotted  against 
the  potential  of  the  collector,  a  straight  line  should  be  obtained.  If  S  is 
the  slope  of  this  line,  one  has 


S  = 


4  A2  I2e 
;rkT 
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or 

I/VT"  =  >/jrk/4e  Vs/A  =  0.00822  >fs/ A  amp  cm"2  deg" 1/2  (10) 

The  intercept  of  the  straight  line  on  the  V-axis  is  obtained  from 
Eq.  (9)  by  placing  i  =  0  and  corresponds  to  a  value  of  V  given  by 

Vi=Q  =  -kT/e  =  -T/ 11,  600  volts 


Thus,  the  straight  line  represented  by  Eq.  (9)  crosses  the  V-axis  at  a 
point  T/ll,  600  volts  negative  with  respect  to  the  potential  of  the  space 
around  the  probe. 


From  the  slope  S  one  may  calculate  directly  the  number  of  ions 
(or  electrons)  per  unit  volume.  The  average  velocity  v  of  the  ions, 
given  by  the  kinetic  theory,  is 


v  = 


8k  T 
ir  m 


and  the  number  of  ions,  n,  corresponding  to  a  current  density  I,  per 
second  per  unit  area  passing  a  plane  is  given  by 


n 


ve 


Substituting,  one  obtains  an  expression  for  the  number  density, 

n  I/e  =  4.  03  x  1013(l/>[T)>/m/me 

and  combining  with  Eq.  (10), 

n  =  --1- —  ^S/Ae  =  3.  32  x  10U(>/S/A)  y}m/ me  (11) 

V  2e/m 

if  yJS  is  expressed  in  amp  volt"^2. 

Near  the  space  potential,  the  square  of  the  current  will  no  longer 
follow  the  straight  line  law  since  the  function  f  is  not  accurately  repre¬ 
sented  by  Eq.  (8)  in  this  region. 

In  the  following  sections,  these,  equations  will  be  applied  to  the 
experimental  data  obtained  with  the  small  platinum  probe  in  a  low 
density  nitrogen  flow.  Since  the  cylindrical  probe  used  in  these  studies 
was  aligned  with  its  axis  parallel  to  the  flow,  the  effects  of  flow  velocity 
can  be  ignored  (Ref.  4). 
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If  the  electrons  in  the  charged  flow  field  have  a  Maxwellian  distri¬ 
bution,  the  electron  temperature  may  be  determined  from  Eq.  (1)  with 
the  electrostatic  probe  being  operated  in  the  retarding  potential  region, 
that  is,  where  the  current- voltage  characteristic  gives  an  indication  of 
the  distribution  of  electron  velocities.  Taking  the  natural  logarithm  of 
both  members  of  Eq.  (1),  one  obtains 

in  i  =  const  +  Ve/(kT)  (12) 

Thus,  if  the  logarithm  of  the  electron  current  i  is  plotted  as  a  function 
of  the  probe  potential  V,  a  straight  line  of  slope  e/(kT)  should  be 
obtained.  Thus,  the  electron  temperature  in  degrees  Kelvin  may  be  de¬ 
termined  from 

T  =  5040 /s  °K  (13) 

where  s  is  the  slope  on  semilog  base  10  paper. 


SECTION  IV 

EXPERIMENTAL  PROCEDURE 


4.1  CALIBRATION  OF  THE  SUPERSONIC  NOZZLE 

The  calibration  of  the  nozzle  exit  flow  was  made  with  a  2.  5-cm- 
diam  pitot  tube  in  the  conventional  manner.  For  various  plenum  pres¬ 
sures,  PQ,  and  stagnation  temperatures,  T0,  the  total  mass  flow 
(Fig.  8)  and  pitot  pressures  were  recorded.  The  inviscid  stagnation 
pressure,  Pq,  was  determined  by  iterating  with  the  pitot  viscous  effect 
function  (Fig.  9)  from  Ref.  5.  The  ratio  Pq/P0  which  is  a  unique  func¬ 
tion  of  the  Mach  number  was  then  determined  and  the  free-stream  Mach 
number  obtained.  Using  the  Mach  number,  plenum  pressure,  and  stag¬ 
nation  temperature  the  free-stream  properties  may  be  calculated. 
Figure  10  presents  the  corrected  values  of  Mach  number  for  various 
values  of  plenum  pressure  for  a  stagnation  temperature  of  290°K  and  a 
nozzle  wall  temperature  of  77°K. 

To  determine  the  variation  of  stream  uniformity,  the  pitot  tube 
was  traversed  across  the  flow  field.  Figure  11  presents  the  pitot  pres¬ 
sure  at  various  positions  across  the  stream  in  the  nozzle  exit  plane  for 
various  stagnation  pressures.  The  variation  of  core  diameter  with 
stagnation  pressure  is  presented  in  Fig.  12. 
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4.2  OPERATIONAL  PROCEDURE  FOR  ION  DENSITY  MEASUREMENT 

Ionized  gas  flow  was  generated  at  various  levels  of  plenum  pressure 
and  microwave  power.  Values  for  relative  microwave  power  were  taken 
from  the  percent-of-power  meters  on  each  generator.  The  generators 
were  individually  set  for  the  same  percent-of-power  output.  A  Tesla 
coil  was  used  to  initiate  ionization  in  each  cavity,  and  once  initiated  it 
was  easily  maintained  over  a  wide  range  of  power.  Tuning  of  the  micro- 
wave  cavities  was  accomplished  with  the  aid  of  a  reflected  power  meter. 
Each  cavity  was  adjusted  for  minimum  standing  wave  ratio. 

In  addition  to  ions,  the  microwave  cavities  also  produce  not  only 
other  excited  species  but  also  dissociation  in  the  test  gas.  The  effects 
of  these  products  in  terms  of  ionospheric  simulation  have  not  been  con- 
*  sidered  in  this  report. 

The  platinum,  cylindrical  probe  previously  described  was  used  to 
obtain  the  voltage -current  characteristic  of  the  ionized  gas  in  the  nozzle 
exit  plane.  The  axis  of  the  probe  was  parallel  to  the  nozzle  axial  center- 
line  and  the  probe  was  traversed  across  the  stream  while  data  were 
taken  at  various  stations. 

Voltage  applied  to  the  probe  was  swept  slowly  (0.01  Hz)  to  avoid 
problems  with  the  d-c  meter  used  to  read  current.  Frequency  response 
of  the  meter  amplifier  is  3  db  down  at  1  Hz.  There  was  a  noticeable  lag 
in  meter  output  if  sweep  rates  greater  than  0.  1  Hz  were  attempted. 


SECTION  V 

EXPERIMENTAL  RESULTS 


5.1  ION  DENSITY  PROFILES 

The  ion  density  was  determined  from  Eq.  (11)  since  the  current- 
voltage  characteristic  was  described  by  Eq.  (9).  Thus,  the  flow  con¬ 
ditions  were  such  that  orbital  motions  of  the  electrons  and  ions  were 
significant.  A  typical  example  of  the  current-voltage  characteristics 
at  various  exit  plane  locations  for  a  plenum  pressure  of  120  MHg  is 
shown  in  Fig.  13.  Negative  values  of  length  are  the  distances  measured 
left  of  the  nozzle  centerline,  whereas  positive  values  are  measured  to  the 
right  of  centerline.  From  these  curves,  graphs  of  the  square  of  the  ion 
current  versus  probe  voltage  were  constructed.  A  complete  set  of  these 
plots  for  the  various  flow  conditions  is  presented  in  Figs.  14  through  22. 
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Since  the  current  is  determined  by  orbital  motion,  the  square  of  the 
current  when  plotted  versus  the  probe  voltage  gives  a  line  nearly  straight. 
The  deviation  of  the  points  from  a  straight  line  is  caused  by  the  arrival 
of  electrons.  The  slopes  of  these  lines  were  determined  and  the  ion  num¬ 
ber  density  calculated.  Normalized  ion  density  profiles  for  various  flow 
conditions  and  microwave  cavity  power  are  presented  in  Figs.  23 
through  26.  Centerline  ion  densities  ranged  from  5.  75  x  10®  to 
1.  69  x  10®  ion/cm^  for  Mach  numbers  of  2.  86  to  3.  17,  respectively. 

In  all  cases,  the  ion  profiles  were  narrower  than  the  corresponding 
pitot  pressure  profiles.  Also,  the  centerline  ion  densities  increased  as 
the  microwave  power  decreased,  as  shown  in  Table  II.  One  possible 
explanation  for  this  increase  is  that  as  the  microwave  power  increases, 
the  efficiency  of  ion  generation  may  decrease  because  of  a  change  in  the 
properties  of  the  plasma,  for  example,  the  plasma  impedance.  This 
decrease  in  ion  generation  rate  would  result  in  smaller  exit  plane  ion 
densities,  thus  giving  the  observed  results. 


5.2  ESTIMATION  OF  THE  ELECTRON  TEMPERATURE 

The  electron  temperatures  were  estimated  using  Eq.  (13)  after 
determining  the  values  of  the  slope,  s,  from  a  semilogarithmic  plot  of 
the  current-voltage  characteristics.  Figures  27  through  29  present 
these  graphs  for  various  flow  conditions.  Table  III  presents  the  estimated 
values  of  the  electron  temperature  at  various  probe  positions.  Table  IV 
gives  T  for  various  Mach  numbers,  and  Table  V  presents  T  for  several 
values  of  microwave  power.  An  examination  of  these  tables  reveals  a 
great  deal  of  scatter  in  the  estimation  of  the  electron  temperature.  For 
the  test  conditions  used  in  this  investigation,  the  value  of  T  should  be 
essentially  constant.  The  scatter  is  attributed  primarily  to  the  difficulty 
in  replotting  the  data  from  the  current-voltage  characteristic  onto  the 
semilogarithmic  graph.  An  average  of  the  estimated  values  of  the  elec¬ 
tron  temperatures  gives  3000°K  for  T.  The  smallest  calculated  value 
of  T  was  2180aK  (Table  V),  whereas  the  largest  value  was  3980°K 
(Table  in).  Thus,  the  value  of  electron  temperature  in  this  report  is 
taken  to  be  3000°K  ±  900°K. 


5.3  CONCLUSIONS 

The  feasibility  of  the  use  of  microwave  cavities  for  the  production  of 
low  density  ionized  flows  has  been  demonstrated.  Measurements  of  the 
characteristics  of  the  neutral  flow  were  carried  out  using  a  pitot  probe, 
and  having  applied  the  appropriate  viscous  corrections,  the  neutral  flow 
Mach  number  and  density  were  determined.  A  cylindrical,  platinum 
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electrostatic  probe  with  its  axis  aligned  parallel  to  the  nozzle  axial 
centerline  was  used  to  record  the  current -volt age  characteristics  of 
the  partially  ionized  gas.  The  flow  conditions  were  such  that  the 
orbital  motions  of  the  ions  and  electrons  were  significant,  and  the 
appropriate  equations  were  derived  and  used  to  reduce  the  experimen¬ 
tal  data.  Centerline  exit  plane  ion  densities  of  5.  75  x  10^  ions/ cm  ^ 
to  1.  69  x  10®  ions/cm^  for  Mach  numbers  of  2.  86  to  3.  17,  respectively, 
were  measured  with  estimated  electron  temperatures  of  3000°K  ±  900°K. 
In  all  cases  investigated,  the  ion  density  profiles  were  narrower  than 
the  corresponding  pitot  pressure  profiles.  Table  VI  presents  a  compar¬ 
ison  of  the  desired  ionosphere  flow  characteristics  and  those  actually 
measured  in  the  experiment.  The  simulation  of  the  ionosphere  flow 
properties  is  good,  as  shown  in  Table  VI,  although  the  experimental 
value  of  the  neutral  density  and  electron  temperature  is  somewhat 
higher  than  those  desired.  This  experimental  study  has  only  been 
preliminary,  and  further  work  is  required  (for  example,  calibration 
using  other  electrostatic  probes  and  sophistication  and  improvement  of 
the  data  acquisition  techniques)  in  order  to  provide  an  ionospheric 
simulation  facility  for  the  ground  testing  of  rocket -borne  atmospheric 
probes,  gages,  and  other  devices. 
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Fig.  1  Typical  Electron  Distributions  in  the  Ionosphere 


15 


os 


Fig.  2  Research  Vacuum  Chamber 
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Fig.  4  Gas  Addition  System  and  Microwave  Generators 
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Fig.  5  Schematic  of  the  Electrostatic  Probe  Instrumentation 
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Distance  from  Nozzle  Centerline,  in. 

Fig.  11  Pitot  Pressure  Profiles  in  Nozzle  Exit  Plane 
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Plenum  Pressure,  MHg 


Fig.  12  Nozzle  Exit  Plane  Core  Diameter  for  Various  Plenum  Pressures 
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Fig.  13  Voltage-Current  Characteristics  for  Various  Probe  Positions 
in  Nozzle  Exit  Plane 
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Fig.  14  l2-V  Curves  for  Mach  Number  3.17  and  100-percent  Microwave  Power 
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Fig.  15  l2-V  Curves  for  Mach  Number  3.17  and  75-percent  Microwave  Power 
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Fig.  17  l2-V  Curves  for  Mach  Number  3.06  and  100-percent  Microwave  Power 
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Fig.  18  l2-V  Curves  for  Mach  Number  3.06  and  50-percent  Microwave  Power 
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Fig.  19  l2-V  Curves  for  Mach  Number  2.91  and  100-percent  Microwave  Power 
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Fig.  20  l2-V  Curves  for  Mach  Number  2.91  and  50-percent  Microwave  Power 
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Fig,  21  l2-V  Curves  for  Mach  Number  2.86  and  100-percent  Microwave  Power 


35 


amp 


AE  DC-TR-71  -254 


0  1 _ L— J _ L _ I _ I _ 1 _ I _ I _ I _ I  1  Y  T  I  I _ 1  I  I  I  I  I _ I  I  I 

-15  -10  -5  0 


V,  volts 

Fig.  22  l2-V  Curves  for  Mach  Number  2.86  and  50-percent  Microwave  Power 
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Fig.  23  Normalized  Exit  Plane  Ion  Density  Profiles  for  IViach  Number  3.17 
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Fig.  24  Normalized  Exit  Plane  Ion  Density  Profiles  for  Mach  Number  3.06 
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Fig.  25  Normalized  Exit  Plane  Ion  Density  Profiles  for  Mach  Number  2.91 
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Fig.  26  Normalized  Exit  Plane  Ion  Density  Profiles  for  Mach  Number  2.86 
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Fig.  27  Nozzle  Exit  Plane  l-V  Curves  for  Various  Probe  Positions 
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Fig.  28  Nozzle  Exit  Plane  l-V  Curves  for  Various  Mach  Numbers 
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Fig.  29  Nozzle  Exit  Plane  1*V  Curves  for  Various  Microwave  Powers 
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TABLE  I 

CHARACTERISTICS  OF  THE  IONOSPHERE 


Height, 

km 

Neutral 

Density 

(cm~3) 

Electron 

Density 

(cm~3) 

Electron 

Temperature, 

°K 

Mean 

Molecular 

Weight 

100 

3  x  1013 

(2-100)  x  103 

230 

28 

300 

3  x  109 

(1-200)  x  105 

1000 

24 

500 

5  x  107 

(4-10)  x  105 

1800 

20 

700 

6  x  106 

(2  -  5)  x  105 

2000 

16 

1000 

105 

105 

3000 

16 

TABLE  li 

CENTERLINE  ION  DENSITIES  IN  THE  NOZZLE  EXIT  PLANE  FOR 
VARIOUS  MICROWAVE  POWERS 


Mach  Microwave  Power, 

Number  percent 

Centerline  Ion  Densities 
(cm"3) 

3.  17 

100 

1.  37 

X 

106 

3.  17 

75 

1.39 

X 

106 

3.  17 

50 

1.  69 

X 

106 

3.06 

100 

1.  06 

X 

106 

3.06 

50 

1.36 

X 

106 

2.  91 

100 

8.  98 

X 

105 

2.  91 

50 

1.02 

X 

106 

2.  86 

100 

5.  75 

X 

105 

2.  86 

50 

7.  09 

X 

105 
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TABLE  III 


ESTIMATED  ELECTRON  TEMPERATURES  IN  THE 
NOZZLE  EXIT  PLANE  FOR  MACH  NUMBER  *  3.17 

Position, 

Temperature. 

in. 

°K 

-7 

3980 

-6 

2960 

-1 

3200 

1 

2760 

3 

3220 

100-percent  Microwave  Power 


TABLE  IV 

ESTIMATED  CENTERLINE  ELECTRON  TEMPERATURES  FOR 
VARIOUS  MACH  NUMBERS  IN  THE  NOZZLE  EXIT  PLANE 


Mach  Temperature, 


Number 

°K 

3. 17 

3200 

3.06 

3700 

2.  91 

2700 

100-percent  Microwave  Power 
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TABLE  V 

ESTIMATED  CENTERLINE  ELECTRON  TEMPERATURES  FOR 
VARIOUS  VALUES  OF  MICROWAVE  POWER  IN  THE  NOZZLE  EXIT  PLANE 


Microwave  Power, 
_ percent _ 


Temperature, 

°K 


100 

75 

50 


3200 

2180 

2800 


Mach  Number  =  3. 17 


TABLE  VI 

COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  DESIRED 
IONOSPHERE  CONDITIONS 


Desired  Conditions  Experimental  Results 


Mach  Numbers  2-5 

Ion  Densities  103  -  107  ions/cm^ 

Electron  Temperature  230  -  1800°K 
Neutral  Density  10^  -  10®  cm"® 


2.  86  -  3.  17 

5.  75  x  105  -  1.69  x  10®  ions/cm3 
3000°K  ±  900°K 
(1.49-2.04)  x  1014  cm-3 
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